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Abstract 


The Tazy or Kazakh National sighthound has been officially recognized as the national heri- 
tage of Kazakhstan. Comprehensive genetic studies of genetic diversity and population 
structure that could be used for selection and conservation of this unique dog breed have 
not been conducted so far. The aim of this study was to determine the genetic structure of 
the Tazy using microsatellite and SNP markers and to place the breed in the context of the 
world sighthound breeds. Our results showed that all 19 microsatellite loci examined were 
polymorphic. The observed number of alleles in the Tazy population varied from 6 (INU030 
locus) to 12 (AHT137, REN169D01, AHTh260, AHT121, and FH2054 loci) with a mean of 
9.778 alleles per locus. The mean number of effective alleles was 4.869 and ranged from 
3.349 f to 4.841. All markers were highly informative (PIC values greater than 0.5) and ran- 
ged from 0.543 (REN247M23 locus) to 0.865 (AHT121 locus). The observed and expected 
heterozygosities in a total population were 0.748 and 0.769 and ranged from 0.746 to 0.750 
and 0.656 to 0.769, respectively. Overall, the results confirmed that the Tazy breed has a 
high level of genetic diversity, no significant inbreeding, and a specific genetic structure. 
Three gene pools underlie the genetic diversity of the Tazy breed. SNP analysis using the 
CanineHD SNP array, which contains more than 170,000 SNP markers, showed that the 
Tazy breed is distinct from other sighthound breeds and genetically related to ancient east- 
ern sighthound breeds sharing the same branch with the Afghan Hound and the Saluki. The 
results, together with archeological findings, confirm the ancient origin of the breed. The find- 
ings can be used for the conservation and international registration of the Tazy dog breed. 


Introduction 


The Tazy is a breed of sighthound common in Kazakhstan. Kazakhs hunted hares, wild boars, 
foxes, badgers, deer, wolves and saigas with the Tazy since ancient times. The role of the Tazy 
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is well documented by famous travelers and explorers of the 18th century such as Pallas P. S., 
Falk I. P., and Georgi I. G. The description of the Tazy as a distinct breed was made by Mache- 
varianov P. M. [1], Bogdanov M. [2], Sabaneev L. P. [3], and Sludsky A.A. [4]. Unfortunately, 
the popularity of the Tazy declined when more modern hunting methods became available. 
Their population size seems to have gradually decreased in recent decades. Unofficial estimates 
indicate that there are currently only about 3,000 Tazy dogs left, but no more than 350 high- 
quality dogs with pedigrees up to the fourth generation. After the resumption of efforts to pre- 
serve this unique breed, the Tazy was granted the status of a national dog breed, and the mod- 
ern standard of the Tazy breed was developed (S1 Table). 

Nevertheless, this unique breed is not recognized by any of the major kennel clubs, neither 
the American Kennel Club (AKC), the United Kennel Club (UKC), the Canadian Kennel 
Club (CKC), nor the Fédération Cynologique Internationale (FCI) and for most dog breeders 
in the world, the Tazy is still unknown and confused with the Saluki. In fact, the phenotypic 
differences between the Tazy and the Saluki are barely perceptible and concern the structure 
of the coat, the shape of the head and the chest. There is a hypothesis that these breeds and 
other eastern sighthounds (Sloughi, Azawakh, Afghan Hound, etc.) are descended from a sin- 
gle common ancestor such as the ancient Egyptian Tezem (or Tesem) [3, 5]. However, 
despite the popularity of this hypothesis, some authors deny that there is evidence that the 
sighthounds are descended from the Tezem [6]. The Sloughi has been cited as a second likely 
common ancestor [5], but mitochondrial DNA analysis revealed that the Tazy has no 
sequences in common with the Sloughi but does have sequences in common with the Saluki 
[7]. This agrees well with the hypothesis that the Tazy probably originated from a cross 
between the Saluki and the Mongolian Shepherd [5], which is supported by breed cluster 
data based on SNP markers showing that the Saluki is one of the oldest breeds in the world 
[8-10]. An alternative hypothesis is that the ancestors of sighthounds were the Tazy or Hor- 
taya Borzaya, since Kazakhstan and the Don and Dnieper steppes are possible centers for the 
domestication of the horse [11]. It is known that the first sighthounds may be found in those 
ancient cultures where the hunter is a horseman, archer, and nomad, i.e., in cultures where 
the horse was domesticated, and tribes moved freely in open landscapes [5]. The origin of the 
sighthounds was probably closely related to the process of domestication of the horse and the 
development of horsemanship, which makes Kazakhstan a possible place of origin of the 
sighthounds [5]. 

So, the discussion about the origin of the sighthounds and the position of the Tazy among 
the other breeds continues, and it becomes important to study the genetic structure of the 
Tazy and to get a picture of how well the breed self-segregates and relates to similar breeds. 
Microsatellite markers, or Short Tandem Repeat (STR) markers, are a well-known, effective, 
and powerful tool commonly used to study the genetic structure and diversity of dog breeds 
[12-15]. Evaluation and monitoring of genetic variability facilitates selection of the optimal 
mating partner and determination of the degree of inbreeding. Information obtained through 
microsatellite marker studies is particularly useful for native dog breeds. Recently, SNP-based 
analysis has opened new possibilities for genetic analysis of dog breeds. A few studies have 
reported the successful identification of breed structure and evolutionary relationships among 
dog breeds using high-density SNP microarrays [9, 15-17]. To our knowledge, the genetic 
diversity and structure of the Tazy dogs has not been studied. The aim of our study was to pro- 
vide insight into the genetic structure and history of the Tazy using two effective and widely 
used genetic markers such as STR loci and SNPs. To determine the position of the Tazy 
among other sighthounds, we characterized the genetic profile of the Tazy by comparing it to 
recognized sighthound breeds and estimated the genetic differentiation within and between 
these breeds. The results of the study could be a unique tool for developing strategies for 
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effective management and conservation of the Tazy’s genetic resources and for international 
recognition of the breed. 


Material and methods 
Sample collection and DNA extraction 


All experimental procedures were approved by the Ethics Committee of the Institute of 
Human and Animal Physiology, Almaty, Kazakhstan (# 3, September 15th, 2020). 

DNA samples were obtained by collecting buccal swabs and/or blood samples at dog shows, 
special events, and mail donations. Between April 2021 and May 2022, a total of 134 samples 
were collected from six breeds of sighthound: Whippet (n = 5), Russian Borzoi (n = 4), Grey- 
hound (n = 8), Afghan Hound (n = 2), Saluki (n = 1), and Tazy (n = 114). The Tazy dogs were 
collected in three locations: South Kazakhstan (TSK, n = 58), North Kazakhstan (TNK, 

n = 46), and East Kazakhstan (TEK, n = 10). Assessment of compliance with the breed stan- 
dard was performed by the expert of the national breed affiliated association "Kansonar". All 
owners gave informed written consent to use samples from their dogs for genetic studies. 
Detailed pedigree information was requested for all dogs. When pedigree information was 
available, only minimally related animals were included in the sample. For those animals for 
which no pedigree was available, non-relatedness was verified by the handlers. DNA was 
extracted using QiaAmp DNA extraction kits according to the manufacturer’s protocol (Qia- 
gen, Valencia, CA). 


STR genotyping 

A total of 19 microsatellite markers (AHTk211, CXX279, REN169018, INU055, REN54P11, 
INRA21, AHT137, REN169D01, AHTh260, AHTk253, INU005, INU030, amelogenin, FH2848, 
AHT121, FH2054, REN162C04 AHTh171, and REN247M23) recommended by ISAG were 
amplified using a commercial kit Canine ISAG STR Parentage Kit (Thermo Fisher Scientific, 
CA, USA) for a total of 134 DNA samples. Genotyping was performed using the SeqStudio™ 
Genetic Analyzer (Thermo Fisher Scientific, CA, USA). Microsatellite alleles were processed 
and manually verified using GeneMapper™ Software 6 (Thermo Fisher Scientific, CA, USA). 


SNP genotyping 

SNP genotyping using an lumina Infinium CanineHD Genotyping BeadChip (Illumina Inc. 
San Diego, CA) was performed for a total of 39 DNAs of the Tazy with the highest expert 
scores, whose functional health and body characteristics meet the standard for the Tazy breed 
and are able to work and perform functions according to the specific characteristics of the 
Tazy breed (S1 Fig). SentrixBarcode and SentrixPosition on the chip are listed in S2 Table. 
Genotype data have been deposited in an Open Science Framework repository and are avail- 
able at: DOI 10.17605/OSF.IO/S5SHWU. 


Analysis of the STR data 


Allele frequencies of 18 STR loci were used to determine average alleles/locus (Na), average 
effective alleles/locus (Ne), observed heterozygosity (Ho), expected heterozygosity (He), 
inbreeding coefficient F, and principal coordinate analysis (PCoA) using GenAIEX 6.5 [18- 
20]. Polymorphism information content (PIC) and deviations from Hardy-Weinberg equilib- 
rium (HWE) were tested using Cervus 3.0.7 [21]. Bonferroni corrections were applied for 
HWE estimation. Bayesian clustering was performed using Structure v2.3.4 [22, 23]. Simula- 
tions were performed assuming a model of admixture and correlated allele frequencies with 10 
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000 burn-in steps and 100 000 MCMC steps. Ten independent runs were performed for each 

assumed number of clusters (K), with K varying from 1 to 10. The most likely number of clus- 
ters was determined in CLUMPAK25 (http://clumpak.tau.ac.il/index.html) using the method 
of Evanno et al. [24]. The analysis was performed twice, the first time using the Tazy database 
and a second time using the database of all sighthound breeds. 


Analysis of the SNP data 


The SNP genotype data generated by the iScan system were loaded into the Illumina Geno- 
meStudio program to perform the primary data analysis and generate a final custom report 
(PED and MAP) for downstream analysis. Data obtained from the 39 Tazy were merged with 
publicly available SNP array data of 89 dogs from seven sighthound breeds and 14 Gray 
Wolves downloaded from the Dryad repository (datadryad.org, doi:10.5061/dryad.v9t5h; 
doi:10.5061/dryad.pm7mt). The sample code and corresponding breed are listed in $3 Table. 
In the PLINK 1 (www.cog-genomics.org/plink/1.9/) [25] Input Report, 166,171 SNPs of the 89 
samples from seven breeds and 14 Gray Wolves and 172,115 SNPs of the 39 samples from the 
Tazy dogs were filtered using the following steps: (1) removal of very closely related individuals 
PI_HUT > 0.4; (2) filtering of SNPs that have an exact Hardy-Weinberg equilibrium (—hwe 
0.01); (3) removal of SNPs on the X and Y chromosomes (—not-chr); (4) selection of only 
SNPs with minor allele frequency (—maf) > 0.05; (5) calling rate SNP (—geno) 0.05; (6) 
removal of SNPs with pairwise genotypic associations (r2) > 0.2 within a window of 50 SNPs 
(—indep-pairwise 50 5 0.2). The number of SNPs retained for calculations after the filtering 
process was 40,229 autosomal SNPs. The PCoA of unrelated dogs was performed using PLINK 
1.9 software and visualized in the R package "ggplot2" [26, 27]. 


Results 
Diversity analysis based on the STR dataset 


A total of 176 alleles were detected at the 18 microsatellite loci. Table 1 shows the variability 
parameters of the analyzed loci in a total population of the Tazy and in three subpopulations. 

The highest number of alleles was detected for AHT137, REN169D01, AHTh260, AHT121, 
and FH2054 (12 alleles for each locus), and the lowest for INU030 (6 alleles per locus). The 
average number of alleles per locus was 9.778 and ranged from 5.056 for TEK to 8.833 for 
TSK. The average number of effective alleles was 4.869 and ranged from 3.349 for TEK to 
4.841 for TSK. There were no significant differences in Ho among all three subpopulations 
(Ho = 0.746-0.750). In contrast, He was highest in TSK (0.769 vs. 0. 0.656 in TEK and 0.759 in 
TNK). The inbreeding coefficient within subpopulations (F) ranged from -0.121 (TEK) to 
0.030 (TSK), whereas the overall inbreeding coefficient for all Tazy dogs was 0.030. The mean 
value of PIC for all Tazy dogs was estimated to be approximately 74%, ranging from 54% for 
REN247M23 to 87% for AHT121 (Table 2). PIC above 60% was observed for all STR markers 
except REN247M23, whereas a value above 80% was observed for 5 loci (AHT137, 
REN169D01, AHT121, FH2054, and AHTh171). 

The Hardy-Weinberg test performed separately for each locus showed no deviation from 
the expected frequencies after Bonferroni correction (P > 0.05), except for INU030 (P < 0.01) 
(Table 3). 


Tazy population structure based on the STR dataset 


To further characterize the genetic variation of breed the STRUCTURE analysis was applied 
(Fig 1). The Delta K results showed that the optimal number of genetic clusters representing 
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Table 2. PIC of 18 STR markers of the Tazy. 
Locus PIC 
AHTk211 0.614 
CXX0279 0.746 
REN169018 0.789 
INU055 0.653 
RENS54P11 0.789 
INRA21 0.702 
AHT137 0.819 
REN169D01 0.861 
AHTh260 0.701 
AHTk253 0.609 
INU005 0.764 
INU030 0.653 
FH2848 0.774 
AHT121 0.865 
FH2054 0.811 
REN162C04 0.765 
AHTh171 0.842 
REN247M23 0.543 
Mean 0.739 
SE 0.094 
Abbreviations: PIC = Polymorphism information content. 
https://doi.org/10.1371/journal.pone.0282041 1002 
Table 3. HWE for 18 STR markers in total population. 
Locus P (with Bonferroni correction) Significance 
AHTk211 0.320 NS 
CXX0279 0.888 NS 
REN169018 0.868 NS 
INU055 0.971 NS 
RENS54P11 0.952 NS 
INRA21 0.083 NS 
AHT137 0.694 NS 
REN169D01 0.836 NS 
AHTh260 0.801 NS 
AHTk253 0.090 NS 
INU005 0.360 NS 
INU030 0.007 s** 
FH2848 0.457 NS 
AHT121 0.903 NS 
FH2054 0.101 NS 
REN162C04 0.811 NS 
AHTh171 0.229 NS 
REN247M23 0.138 NS 
Abbreviations: NS = Not significant; 
** P <0.01, 
https://doi.org/10.1371/journal.pone.0282041 .t003 
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TSK TNK TEK 


Fig 1. Genetic structure of the Tazy subpopulations. Bayesian clustering on the STR dataset of 114 dogs performed with STRUCTURE v2.3.4 after 
correction by Evanno et al. (CLUMPAK)): A- Delta K results; B and C- bar plots where each dog accession is represented by a single vertical line, and 
this line shows colored segments representing the relative percentage of membership in the K cluster. B—bar plot grouped by clusters; C—bar plot 


grouped by subpopulations. Abbreviations of subpopulations: TSK—Tazy of South Kazakhstan, TNK—Tazy of Northern Kazakhstan, TEK—Tazy of 
Eastern Kazakhstan. 


https://doi.org/10.1371/journal.pone.0282041 .g001 


the most similar individuals was K = 3 (Fig 1A), indicating that three gene pools shaped the 
genetic structure of the Tazy population under study (Fig 1B). Based on K = 3, the proportion 
of each of the three gene clusters is present in each subpopulation (Fig 1C). S4 Table shows the 
detailed results on the composition of the clusters and the percentage of membership in each 
accession for three clusters in one run. The distance between each cluster is shown in Table 4. 
Cluster 3 is very similar to clusters 1 and 2 (0.040 for both), while clusters 1 and 2 differ slightly 
more (0.082). The average distances (expected heterozygosity) between individuals in the same 
cluster were found to be highest for cluster 3 (0.808) (Table 4). The estimated mean Fst value 
for accessions in cluster 1 was 0.164, in cluster 2-0.157, and in cluster 3-0.003. 

The PCoA of the STR data across three subpopulations revealed three significant axes: axis 
one explained 5.58% of the variation present, axis two 4.93%, and axis three 4.65%. On both 
axis one and axis two (Fig 2), the TSK and TNK subpopulations had the greatest variation in 


allelic diversity. All three subpopulations were mixed, confirming the results of Bayesian 
clustering. 


Breed relationships based on the STR dataset 


Then we applied STRUCTURE to STR data from the Tazy and other five sighthound breeds to 
analyze the relationships. Plots for K values from 1 to 10 are shown in Fig 3A. The software 
identified K = 5 as the most likely number of genetic clusters in our samples (Fig 3B). $5 Table 
shows the detailed results on the composition of the clusters and the percentage of 


Table 4. Allele-freq. divergence (Net nucleotide distance), computed using point estimates of P through Structure 


analysis. 

Cluster 1 2 3 Average distances Mean value of Fst 
1 ~ 0.082 0.040 0.689 0.164 

2 0.082 - 0.040 0.715 0.157 

3 0.040 0.040 - 0.808 0.003 


Abbreviations: Fst = Degree of differentiation among populations 


https://doi.org/10.1371/journal.pone.0282041 .t004 
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Fig 2. PCoA plot of 114 Tazy dogs from three different subpopulations based on STR data. 
https://doi.org/10.1371/journal.pone.0282041 .g002 


membership in each accession for five clusters in one run. In the visualized structure of one 
run at K = 5, the Greyhound, the Afghan Hound, the Whippet, and the Russian Borzoi form 
distinct clusters consisting of dogs of the respective breeds (Fig 3A). At K = 5, the contribution 
of the Tazy and the Saluki to the Afghan Hound can be detected, but at K = 6, the Afghan 
Hound separates and remains relatively stable. The Tazy and the Saluki show similar mixed 
profiles (Fig 3A). Similarly, in the Structure output presented as a triangle plot (Fig 3C), the 
Tazy and the Saluki are close together in one corner and are assigned to one cluster. The dis- 
tance between each cluster is shown in Table 5. The largest distance was found for cluster 2 
Greyhound vs cluster 5 Russian Borzoi, (0.216), followed by cluster 2 vs cluster 3 Afghan 
Hound (0.197). Cluster 1 Tazy and Saluki was related to clusters 3 and 4 Whippet (0.082 and 
0.085, respectively), as shown in the tree plot (Fig 3D). The mean Fst value (Table 5) for the 
accessions in cluster 1 Tazy and Saluki had the lowest estimated value compared to other 
breeds (0.003). 

The PCoA (Fig 4), both the first and second axes, which accounted for 11.74% and 9.16% of 
the total variance, respectively, separated the Saluki from the Tazy dogs. The Tazy and the 
Afghan Hound were clustered together, which is consistent with the result of the STRUC- 
TURE analysis, in which the relatedness of these two breeds were indicated. 


Breed relationships based on the SNP dataset 


Since Bayesian clustering of the STR dataset did not reveal any divergence between the Tazy 
and the Saluki, we performed SNP-based analyses for 39 Tazy dogs. The resulting dataset was 
merged with publicly available SNP data for 89 dogs from seven sighthound breeds and 14 
Grey Wolves. Relationships between breeds were visualized using PCoA and the phylogenetic 
dendrogram. The PCoA plot shows stratification between breeds and no significant hybridiza- 
tion of sighthounds with wolves (Fig 5). The eight sighthounds formed eight groups. There 
was a clear genetic split between Tazy and Saluki. The Afghan Hound was found between the 
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Tazy and the Saluki, confirming the results of the tree plot and PCoA based on the data from 
STR. The Russian Borzoi, the Bloodhound, the Otterhound, the Whippet and the Greyhound 
were closer and had a distance to the Tazy, the Afghan Hound and the Saluki. 


In the phylogenetic tree, the Tazy comes from the same branchpoint with the Afghan 
Hound and the Saluki (Fig 6), showing that there is an ancient common ancestry between 


Table 5. Allele-freq. divergence (Net nucleotide distance), computed using point estimates of P through Structure 


analysis. 

Cluster 1 3 4 Average distances Mean value of Fst 
1 = 0.107 0.082 0.085 0.110 0.812 0.003 

2 0.107 = 0.197 0.192 0.216 0.609 0.280 

3 0.082 0.197 - 0.142 0.193 0.628 0.337 

4 0.085 0.192 0.142 - 0.195 0.623 0.305 

) 0.110 0.216 0.193 0.195 - 0.575 0.354 


https://doi.org/10.1371/journal.pone.0282041 .t005 
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these three breeds, but a divergence between the Tazy/Afghan Hound and the Saluki, and then 
a more recent divergence between the Tazy and the Afghan Hound. Similarly, the Greyhounds 
and the Whippets share a common ancestral gene pool, as well as the Otterhound and the 
Bloodhound. 


Discussion 


This is the first study describing the genetic structure of the Tazy breed. We focused on three 
regions of Kazakhstan, where the Tazy has been an important hunting dog for many centuries, 
to evaluate the parameters of genetic diversity. Using the 18 STR-marker panel, more than 100 
Tazy dogs were genotyped. The results show that the breed has a high level of genetic diversity. 
The number of alleles that have the greatest genetic influence on heterozygosity, known as Ne, 
averaged 4.869 in the general population, which is higher than the Ne reported for other native 
dog breeds: 3.0-3.5 for the Italian Greyhound [28], 3.3 for the German Shepherd, 3.5 for the 
Maltese, 3.2 for the Biewer Yorkshire Terrier, 3.5 for the Yorkshire Terrier [29]. All markers 
were highly informative (PIC values greater than 0.5) [30], suggesting the efficacy of these 
markers for the assessment of genetic variation among the Tazy dogs. The average polymor- 
phism value based on all Tazy dogs was 0.739, a value very similar to that obtained for Grey- 
hounds from South Korea (0.73) [13]. In contrast, the PIC value in a Polish Greyhound 
population was much lower (0.597), as reported by Goleman et al. [31]. The highest level of 
polymorphism was observed at the REN169D01 and AHT121 loci, where the PIC value 
exceeded 0.850. An equally high degree of polymorphism of these markers was found in 28 
Borzoi dogs in Poland (PIC = 0.756 and 0.732, respectively) [32]. The lowest level of polymor- 
phism was observed for the marker REN247M23, with a PIC value of 0.543. A study of 28 bor- 
zoi dogs also revealed one of the lowest levels of polymorphisms of this marker, with a PIC 
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39 Tazy dogs, 74 dogs from seven sighthound breeds and 14 Grey Wolves. 
https://doi.org/10.1371/journal.pone.0282041.g005 


value of 0.590 [32]. We also found the lowest heterozygosity for this marker (Ho = 0.518). The 
mean observed heterozygosity calculated for all STR was 0.748. A similar degree of heterozy- 
gosity was found for Jack Russell Terriers (Ho = 0.75), Yorkshire Terriers (Ho = 0.73), and 
mixed breeds (Ho = 0.73) [33]. A study of Italian Greyhounds found lower Ho values in the 
range of 0.60-0.62 [28]. The estimated He values were slightly higher than the Ho values for 
the whole population and for the Tazy of the Southern and Northern region. Thus, the fixation 
index values were positive, indicating inbreeding, while the mean F value for the Tazy of the 
Eastern region had a negative value (-0.121), indicating that the dogs were not inbred. How- 
ever, the detected inbreeding was low and did not exceed the value considered extremely high 
(> 0.1) [34]. 

The obtained high degree of diversity and allelic richness suggests that the Tazy breed 
evolved from a relatively large and diverse founder population. Moreover, the intra-breed sub- 
structure of the Tazy identified in the STRUCTURE bar plot may be the result of different 
ancestors. The high diversity inherited from the ancestors was maintained despite possible 
genetic bottlenecks in the last century, when the Tazy lost importance for hunting, and may be 
explained by gene flows between populations. The STRUCTURE bar plot clearly shows migra- 
tion events and exchange of genetic material between all three regions, which can be explained 
by the historically evolved lifestyle of the Kazakhs. The semi-nomadic way of life with seasonal 
migrations of people with their dogs favored a dynamic gene flow, which was partially 
restricted by landscape barriers. However, despite high genetic diversity, the Tazy showed low 
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differentiation when the Tazy dogs were compared with purebred sighthounds (Fst = 0.0033), 
indicating weak effects of reproductive isolation. Similarly, low genetic differentiation and 


high diversity were also found for Portuguese native breeds and were interpreted by the fact 
that these breeds are not closed populations [35]. 


Of particular interest was the comparison between the Tazy and the phenotypically similar 
Saluki. Bayesian clustering on a microsatellite data set showed that by the third cluster 


PLOS ONE | https://doi.org/10.1371/journal.pone.0282041 


March 8, 2023 12/16 


PLOS ONE 


The Tazy dog breed 


calculation, the Tazy breed already differs from the Greyhound, Afghan Hound, Whippet, and 
Russian Borzoi, but has an identical profile to the Saluki. Parker concludes that breeds that do 
not form distinct clusters by microsatellite analysis, but tend to cluster with another breed, 
have only recently been separated [9]. We then used the HD SNP array to create a high-resolu- 
tion population structure of the Tazy and compared the obtained data in a broader context of 
worldwide sighthound breeds. The Tazy appeared to be genetically distant from all other 
breeds included in our study and was part of the eastern sighthound group along with the 
Saluki and Afghan Hound. There was a clear separation between the Tazy and the Saluki dogs. 
In the phylogenetic dendrogram, the Tazy, the Afghan Hound and the Saluki breeds shared 
one clade, and the Tazy formed a subclade with the Afghan Hound. The close relationship of 
the Tazy to the Afghan Hound was to be expected given the phenotypic similarity. A local 
name for the Afghan Hound, Sag-e-Tazi, suggests a possible common ancestry with the Tazy. 
Interestingly, in the phylogenetic dendrogram, only the Saluki was previously closely related to 
the Afghan Hound [8, 9, 17]. These two breeds are considered to be an "ancient" breed based 
on several genetic studies [9, 36] and are believed to have originated > 500 years ago [17, 37]. 
Apparently, the Tazy is the missing link in this group of ancient breeds. The hypothesis that 
the Tazy is a dog with a centuries-old history is consistent with archaeological findings. 
Numerous images of Tazy-like dogs on petroglyphs in South Kazakhstan date from different 
historical periods, up to the 10th-12th centuries [5]. 

Moreover, the results of other significant branching were in clear agreement with the 
known history of the designated breed. In our dendrogram, the Otterhound and the Blood- 
hound share a common ancestral gene pool. The Otterhound is an ancient breed that evolved 
in England from Bloodhounds or other southern hunting dogs, which could be the common 
ancestor of these two breeds [38]. Other breeds that were grouped together are the Greyhound 
and the Whippet, which is consistent with the Parker results [9]. It is believed that the breed 
evolved from terriers and small English Greyhounds and was later crossed with Italian Grey- 
hounds to give the Whippet a sleek appearance [39]. 

The strong limitation of the study is that the number of samples obtained from sighthound 
breeds for the STR analysis was very small due to the lower popularity of these breeds in our 
country. We are aware that a small number of samples should not be used to achieve adequate 
power in admixture mapping studies. With a smaller sample size, a study may not detect small 
or moderate effects [40]. However, in this study, the STR analysis was able to distinguish the 
Greyhound, Afghan Hound, Whippet, and Russian Borzoi with a small number of samples. 
This is evidence that genetic clustering algorithms implemented in STRUCTURE can be suc- 
cessfully used to characterize individuals and populations previously separated with minimal 
genetic contact over hundreds of years. A few studies support the successful clustering of 
breeds based on only a few representatives of each breed [9, 41, 42]. In addition, we present 
STRUCTURE and PCoA results based on STR markers alongside PCoA results based on data- 
sets of 170,000 marker genotyping arrays to test whether the assumptions of the model are 
likely to hold and to validate certain features of the results. We hope that this study is the 
beginning of a detailed demographic and historical analysis of the Tazy dog breed. Future 
studies on the Tazy should focus on the comprehensive study of mtDNA, genetic selection of 
hunting success and determination of disease genes. 


Conclusion 


The Tazy breed consists of individuals with specific morphological and behavioural character- 
istics and represents an undeniable genetic and cultural heritage. The history of the Tazy and 
their role as hunters are well documented. In this study, we used data from STR and SNP 
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sources to study the genome structure of the Tazy breed. We characterized the Tazy as a 
unique breed with a specific population structure. Genetic analysis shows that the Tazy have 
high genetic diversity. One of the most important results of our genetic study is the identifica- 
tion of the Tazy as a genetically divergent ancient dog breed with a strong position in the phy- 
logenetic tree. These results provide the first scientific basis for genetic data to improve and 
maintain this breed in the future and support the proposal for international recognition of the 
Tazy. 


Supporting information 


S1 Fig. Pictures of studied Tazy. 
(PDF) 


$1 Table. Breed standard of the Tazy breed (the Republican Federation of Public Associa- 
tions of Hunters and Hunting Societies "Kansonar" and the Ministry of Agriculture of the 
Republic of Kazakhstan). 

(PDF) 


$2 Table. SentrixBarcode and SentrixPosition of the Tazy samples on the Illumina Infi- 
nium CanineHD Genotyping BeadChip. 
(PDF) 


$3 Table. The sample code and corresponding breed (datadryad.org). 
(PDF) 


S4 Table. Tazy dogs and cluster assignment. 
(PDF) 


$5 Table. Sighthound dogs and cluster assignment. 
(PDF) 


Acknowledgments 


We would like to thank the dog breeders and owners who provided us with samples and infor- 
mation about this unique breed. Our special thanks to the veterinarians and members of "Kan- 
sonar" who helped with sampling and provided expert guidance, and to Dossybayev Kairat for 
his critical contributions. 


Author Contributions 
Conceptualization: Anastassiya Perfilyeva, Kira Bespalova. 


Data curation: Anastassiya Perfilyeva, Kira Bespalova, Sergey Bespalov, Yelena Kuzovleva, 
Zhassulan Zhaniyazov, Olga Vishnyakova, Inna Nazarenko, Bakhytzhan Bekmanov. 


Formal analysis: Anastassiya Perfilyeva, Yelena Kuzovleva, Zhassulan Zhaniyazov. 
Investigation: Kira Bespalova, Mamura Begmanova, Yelena Kuzovleva. 


Methodology: Anastassiya Perfilyeva, Kira Bespalova, Mamura Begmanova, Yelena Kuzov- 
leva, Bakhytzhan Bekmanov. 


Project administration: Anastassiya Perfilyeva. 
Resources: Kira Bespalova, Sergey Bespalov, Olga Vishnyakova, Inna Nazarenko. 


Software: Zhassulan Zhaniyazov, Bakhytzhan Bekmanov. 


PLOS ONE | https://doi.org/10.1371/journal.pone.0282041 March 8, 2023 14/16 


PLOS ONE 


The Tazy dog breed 


Supervision: Anastassiya Perfilyeva. 


Visualization: Kira Bespalova, Sergey Bespalov, Bakhytzhan Bekmanov. 


Writing - original draft: Anastassiya Perfilyeva. 


Writing - review & editing: Kira Bespalova, Yuliya Perfilyeva, Ozada Khamdiyeva, Bakhytz- 


han Bekmanov. 


References 


1. 


ga Pen 


13. 


14. 


15. 


16. 


17. 


18. 
19. 


20. 


21. 


22. 


Mayesapnahos [1]. 3anuckn ncoBaro OxOTHUKa Cum6upckon ry6epHun. Mockea: Tun. A. Topjeukaro u 
M. Tepexosa; 1876. 


Borgavos M. Ta3biii uv Kuprickaa 6op3ble. Mpupoga u Oxota 1878;2. 
Ca6anees JI. Co6akuv oxoTHUUbM. Bop3bie u roHuve. Mocksa: Teppa; 1987. 
Cayackui A. Ka3axckaa 60paa Ta3bI u OxoTa c Held. Anma-Ata: KaitHap; 1965. 


Taxos KH, Maxosa AC. ApesHaa uv copeaHeBekoBaa UCTOPUA NOPOAb! GOp3bIx CObaK - Ka3aXCKaA Ta3bl: 
OnbIT peKOHcTpyKuNn. BectHuk PIA3Y 2017; 23:15-25. 


Hawkins R. SIGHTHOUND IDENTITY. The Performance Sighthound Journal USA 2007; 4:36—47. 


Crapon de Caprona D, Fritzsch B. Sloughi, Saluki, Saluqi... Genetic Data Help Separate Semantics 
From Evidence. Hound Breeds Issue 2004. 


Vonholdt BM, Pollinger JP, Lohmueller KE, Han E, Parker HG, Quignon P, et al. Genome-wide SNP 
and haplotype analyses reveal a rich history underlying dog domestication. Nature 2010; 464. https:// 
doi.org/10.1038/nature08837 PMID: 20237475 


Parker HG. Genomic analyses of modern dog breeds. Mammalian Genome 2012; 23. hitps://doi.org/ 
10.1007/s00335-01 1-9387-6 PMID: 22231497 


Hole F, Wyllie C. The Oldest Depictions of Canines and a Possible Early Breed of Dog in Iran. Paléori- 
ent 2007; 33. https://doi.org/10.3406/paleo.2007.5213 


Muxanbckaan A. MATEPUAJIbI O BOP3bIX 2017. 


Bigi D, Marelli SP, Randi E, Polli M. Genetic characterization of four native Italian shepherd dog breeds 
and analysis of their relationship to cosmopolitan dog breeds using microsatellite markers. Animal 
2015; 9. https://doi.org/10.1017/S1751731 115001561 PMID: 26245492 


Kang BT, Kim KS, Min MS, Chae YJ, Kang JW, Yoon J, et al. Microsatellite loci analysis for the genetic 
variability and the parentage test of five dog breeds in South Korea. Genes Genet Syst 2009; 84:245- 
51. https://doi.org/10.1266/ggs.84.245 PMID: 19745574 


Vychodilova L, Necesankova M, Albrechtova K, Hlavac J, Modry D, Janova E, et al. Genetic diversity 
and population structure of African village dogs based on microsatellite and immunity-related molecular 
markers. PLoS One 2018; 13. https://doi.org/10.1371/journal.pone.0199506 PMID: 29940023 


Janes M, Zorc M, Cubric-Curik V, Curik 1, Dove P. Population structure and genetic history of Tibetan 
Terriers. Genetics Selection Evolution 2019; 51. https://doi.org/10.1186/s1271 1-019-0520-4 PMID: 
31881816 


Gajaweera C, Kang JM, Lee DH, Lee SH, Kim YK, Wijayananda HI, et al. Genetic diversity and popula- 
tion structure of the Sapsaree, a native Korean dog breed. BMC Genet 2019; 20. https://doi.org/10. 
1186/s12863-019-0757-5 PMID: 31382890 


Choi L, An JY. Genetic architecture of autism spectrum disorder: Lessons from large-scale genomic 
studies. Neurosci Biobehav Rev 2021; 128. https://doi.org/10.1016/j.neubiorev.2021.06.028 PMID: 
34166716 


Wright S. Coefficients of Inbreeding and Relationship. Am Nat 1922; 56. https://doi.org/10.1086/279872 


Peakall R, Smouse PE. GenALEx 6.5: Genetic analysis in Excel. Population genetic software for teach- 
ing and research-an update. Bioinformatics 2012; 28:2537-9. https://doi.org/10.1093/bioinformatics/ 
bts460 PMID: 22820204 


Peakall R, Smouse PE. GENALEX 6: Genetic analysis in Excel. Population genetic software for teach- 
ing and research. Mol Ecol Notes 2006; 6. https://doi.org/10.1111/j.1471-8286.2005.01155.x 


Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates 
genotyping error increases success in paternity assignment. Mol Ecol 2007; 16:1099—106. https://doi. 
org/10.1111/j.1365-294X.2007.03089.x PMID: 17305863 


Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype 
data. Genetics 2000; 155. https://doi.org/10.1093/genetics/155.2.945 PMID: 10835412 


PLOS ONE | https://doi.org/10.1371/journal.pone.0282041 March 8, 2023 15/16 


PLOS ONE 


The Tazy dog breed 


23. 


24. 


25. 


26. 
27. 


28. 


29. 


30. 


31. 


32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: 
Linked loci and correlated allele frequencies. Genetics 2003; 164. https://doi.org/10.1093/genetics/164. 
4.1567 PMID: 12930761 


Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software 
STRUCTURE: A simulation study. Mol Ecol 2005; 14:261 1-20. https://doi.org/10.1111/j.1365-294Xx. 
2005.02553.x PMID: 15969739 


Chang CC, Chow CC, Tellier LCAM, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: Rising 
to the challenge of larger and richer datasets. Gigascience 2015; 4. https://doi.org/10.1186/s13742- 
015-0047-8 PMID: 25722852 


Wickham H. Package ‘ggplot2‘: Elegant Graphics for Data Analysis. Springer-Verlag New York 2016. 


Galili T. dendextend: An R package for visualizing, adjusting and comparing trees of hierarchical clus- 
tering. Bioinformatics 2015; 31. https://doi.org/10.1093/bioinformatics/btv428 PMID: 26209431 


Pedersen NC, Liu H, Leonard A, Griffioen L. A search for genetic diversity among Italian Greyhounds 
from Continental Europe and the USA and the effect of inbreeding on susceptibility to autoimmune dis- 
ease 2015. https://doi.org/10.1186/s40575-015-0030-9 PMID: 26526059 


Radko A, Podbielska A. Microsatellite dna analysis of genetic diversity and parentage testing in the pop- 
ular dog breeds in Poland. Genes (Basel) 2021; 12. https://doi.org/10.3390/genes12040485 PMID: 
33810589 


Botstein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using 
restriction fragment length polymorphisms. Am J Hum Genet 1980; 32:314. PMID: 6247908 


Goleman M, Balicki |, Radko A, Rozempolska-Ruci nska I, Zi G. Pedigree and Molecular Analyses in 
the Assessment of Genetic Variability of the Polish Greyhound 2021; 11:353. https://doi.org/10.3390/ 
ani11020353 PMID: 33572519 


Radko A, Stota E. Application of 19 microsatellite DNA markers for parentage control in Borzoi dogs. 
Pol J Vet Sci 2009; 12. 


Mellanby RJ, Ogden R, Clements DN, French AT, Gow AG, Powell R, et al. Population structure and 
genetic heterogeneity in popular dog breeds in the UK. Veterinary Journal 2013; 196. htips://doi.org/10. 
1016/j.tvjl.2012.08.009 PMID: 23084740 


Bannasch D, Famula T, Donner J, Anderson H, Honkanen L, Batcher K, et al. The effect of inbreeding, 
body size and morphology on health in dog breeds. Canine Med Genet 2021; 8. https://doi.org/10.1186/ 
s40575-021-00111-4 PMID: 34852838 


Pires AE, Amorim IR, Ginja C, Gomes M, Godinho |, Simées F, et al. Molecular structure in peripheral 
dog breeds: Portuguese native breeds as a case study. Anim Genet 2009; 40. https://doi.org/10.1111/j. 
1365-2052.2009.01849.x PMID: 19298456 


Boyko AR, Boyko RH, Boyko CM, Parker HG, Castelhano M, Corey L, et al. Complex population struc- 
ture in African village dogs and its implications for inferring dog domestication history. Proc Natl Acad 
Sci U S A 2009; 106. https://doi.org/10.1073/pnas.0902129106 PMID: 19666600 


Akey JM, Ruhe AL, Akey DT, Wong AK, Connelly CF, Madeoy J, et al. Tracking footprints of artificial 
selection in the dog genome. Proc Natl Acad Sci U S A 2010; 107. https://doi.org/10.1073/pnas. 
0909918107 PMID: 20080661 


Bell J, Cavanagh K, Tilley L, Smith F. Otterhound. Veterinary Medical Guide to Dog and Cat Breeds, 
2012. https://doi.org/10.1201/b16185-122 


whippet | breed of dog | Britannica n.d. https://www.britannica.com/animal/whippet (accessed Septem- 
ber 30, 2022). 


Lawson DJ, van Dorp L, Falush D. A tutorial on how not to over-interpret STRUCTURE and ADMIX- 
TURE bar plots. Nat Commun 2018; 9. https://doi.org/10.1038/s41467-018-05257-7 PMID: 30108219 


Shannon LM, Boyko RH, Castelhano M, Corey E, Hayward JJ, McLean C, et al. Genetic structure in vil- 
lage dogs reveals a Central Asian domestication origin. Proc Natl Acad Sci U S A 2015; 112. https://doi. 
org/10.1073/pnas.1516215112 PMID: 26483491 


Parker HG, Dreger DL, Rimbault M, Davis BW, Mullen AB, Carpintero-Ramirez G, et al. Genomic Anal- 
yses Reveal the Influence of Geographic Origin, Migration, and Hybridization on Modern Dog Breed 
Development. Cell Rep 2017; 19. https://doi.org/10.1016/j.celrep.2017.03.079 PMID: 28445722 


PLOS ONE | https://doi.org/10.1371/journal.pone.0282041 March 8, 2023 16/16 


